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Abstract 



We obtain the spectrum of glueball masses for the Af = 1 non-conformal cascade 
theory whose supergravity dual was recently constructed by Klebanov and Strassler. The 
glueball masses are calculated by solving the supergravity equations of motion for the 
dilaton and the two-form in the deformed conifold background. 
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1 Introduction 



The original AdS/CFT correspondence |[T|]-[0] was generalized in for branes at con- 

ical singularities (see also |01-|T^ for related work). If instead of locating D3-branes in a 
fiat transversal space they are placed at the vertex of a six dimensional cone with base 
a five dimensional Einstein manifold X5, one is lead to conjecture that type IIB string 
theory on AdS^ x X5 is dual to the low energy limit of the worldvolume theory on the 
D3-branes at the singularity. In particular, a set of D3-branes at the conifold singu- 
larity results on a A/" = 1 superconformal field theory with SU{N) x SU (N) gauge group 
0. Conformal invariance can be broken using fractional D3-branes which are allowed 
to appear in certain singular spaces Fractional D3-branes arise from D5-branes 



wrapped at the collapsed two-cycles of the singularity. In particular, the addition of M 
fractional branes at the singular point of the conifold modifies the gauge group of the 



field theory to SU{N + M) x SU{N). This theory was first investigated in [Tl] to leading 



order in M/N. This solution was then completed to all orders in [18|. However, the 
result found by the authors of [|I^] contains an infrared singularity. In |TP[ , Klebanov and 
Strassler constructed a non singular solution valid from the ultraviolet to the infrared. 
They showed that the theory undergoes a series of Seiberg dualitites as the gauge group is 
sucesively broken. Far in the infrared the D3-branes disappear, and the theory is A/" = 1 
SU{M) with no matter; this theory exhibits confinement, domain walls and screening. 



Using the supergravity description of gauge theories, as in ||2^-[^, we will in the present 
letter find glueball masses for A/" = 1 SU{M) Yang-Mills. In section 2 we will review 
some of the aspects of the solution constructed in |]19| , emphasizing those relevant to our 



calculation. In section 3 we analyze the equations of motion of type IIB supergravity, 
and their final form after dimensional reduction on a sphere. The equations for the type 
IIB dilaton and the two-form fields are numerically solved in section 3 to determine the 
mass spectra of the corresponding glueball modes. Finally, in section 4, we present some 
conclusions and comment on future perspectives^. 



^After completion of this paper, reference |p3| appeared, which partially overlaps with the present 
work. 
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2 Review of Klebanov-Strassler's Solution 



The solution recently found by Klebanov and Strassler arises from the study of D3-branes 
at a singular space. Originally, the authors of ^ studied the conformal field theory on 
D3-branes at a Calabi-Yau singularity dual to type IIB on a AdS^ x T^^ background. 
The corresponding gauge theory is A/" = 1 supersymmetric with SU (N) x SU (N) gauge 
group, and matter content Ai,Bi where i = 1,2. The chiral fields transform as (N, N) and 
(N, N) respectively. The superpotential isW = Ae^-^e^/Tr AiB^AjBi . In the presence of M 
fractional branes the superpotential and matter content are the same but the gauge group 
changes to SU {N + A^) x SU (N) . The chiral superfields are now in the representation 
(N + M, N). The supergravity equations corresponding to this situation were solved, to 
leading order in M/N, in [0], where the relative gauge coupling g^^ — (7^^ was shown 



to run logarithmically. This approximation was completed to all orders in In this 

solution a logarithmic harmonic function warps the conifold. 



,2 

ds 



^ dxndxn + ^ ^^""i-^^^'^ t/r^ + L'^^/\n{r/rs)dslll. (2.1) 



r 



These fractional D3-branes at the singularity are D5-branes wrapped over the collapsed 
S"^ of T^^. D5-branes are sources for the magnetic R-R three-form flux through the S^ 
cycle of T^^ and thus adding M fractional D3-branes implies that the supergravity dual 
of this field theory will involve M units of three-form flux, 

F3 = M, (2.2) 

in addition to units of five-form flux coming from the D3-branes, 

F5 = N. (2.3) 
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This non vanishing three-form is responsible of the conformal symmetry breaking. As 
a consequence of this flux, the two-form B2 is no longer constant; it develops a radial 
dependence [O, 



/ ^2 ~ Me'^ln(r/r,), (2.4) 



while the dilaton remains constant. 
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The supergravity solution obtained in was completed in [|I8| taking into account 



the back reaction of = dB2 and F3 on other fields. This exact solution exhibits an 
effect which was hidden to leading order in M: since F5 = dC^ + B2AF3, F5 also acquires 
a radial dependence, 

F5 = ^5 + *^5, (2.5) 

where 

= /C(r)vol (T^^) = (A^ + agsM^ ln(r/ro))vol (T^^), (2.6) 

with a a constant of order unity. But there is a novelty arising from this solution: the 
five-form fiux present at the ultraviolet scale r = vq may disappear once we reach the 
scale r = r, where IC{f) = 0. This phenomenon can be related to the fact that the fiux 
B2 is not a periodic variable in the supergravity solution, because as this fiux goes 
through a period, /C(r) — > /C(r) — M, so that the five-form fiux is decreased by M units. 
This decrease represents a renormalization group cascade, that was identified in [O] as a 



form of Seiberg duality ||24 



However, the metric ( |2.1| ), representing the logarithmic renormalization group cascade, 
contains a naked singularity at r = r,,, which is the point where the harmonic function 
vanishes, h{rs) = 0. From a physical point of view the singularity represents the end 
point of the cascade as the theory fiows to the infrared, because negative values of are 
unphysical. And it is this singularity in the metric ( p.l| ), that demands a modification, 
at least in the infrared, of the solution. This non singular correction was found in ||T9| , 
where it was argued that the conifold should be replaced by a deformed conifold, 

4 

Y^^i =-'^d^UjZi,=e\ (2.7) 

i=l 

with the singularity removed through the blow up of the in T^^. A simple argument 
in favor of this suggestion comes from the origin of this singularity; it can be related to 
the divergent energy of the F3 field. As there are M units of fiux of F3 through the in 
T^^, when shrinks to zero size F3 diverges. However, if is kept of finite size, as in 
the deformed conifold, there is no need for F3 to diverge. 

A deeper argument comes from a detailed field theory analysis, which shows that the 
spacetime geometry should be modified by the strong dynamics of the infrared limit of 
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the field theory. The authors of showed how the U{1) (Z2M, to be more precise) 
R-symmetry is broken to a Z2 symmetry. This is indeed the symmetry left unbroken on 
the deformed conifold ( p.7| ), Zk — ^ —Zk, instead of the original U{1) Zk — > e'^Zk- The 
metric was then shown to be of the form 



ds 
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h-^'^{T)dXndXn + h^'^{T)dl 



6' 



(2i 



with ds\ the metric of the deformed conifold. In the basis {r, g'^ ^'■■■'^(■?/', 6*1, 62, (pi, (P2)} of 
reference I^Sj this metric becomes diagonal, 



(2.9) 



where 



K(t) 



(sinh(2r) - 2r) 



1/3 



2V3sinh(r) ' 

The harmonic function in (|2.8|) is given by the integral expression 



(2.10) 



h{T) = a 



22/3 



dx (sinh(2x) - 2x)^/^ 



:2.111 



sinh X 

which cannot be evaluated in terms of elementary or special functions. The constant a is 
The solution contains a five-form and three-form fiux. F5 = JF5 + -kj^^ is given by; 



J-5 = gsMH{T)g^ Ag^Ag^Ag^Ag' 



and 



^^•5 = gsM^ 



2l{r) 



15K2(r)/i2(r) sinh2(r)e8/3 
while the three-form is 



A dx^ A dx^ A dx^ A dr. 



G3 — -^3 H 

9s 

= M{g'' A A / + d[F{T){g' A / + / A /)] 
+id[f{r)g'Ag' + k{r)g'Ag']}, 



(2.12) 



(2.13) 



(2.14) 



(2.15) 
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with the functions 



, , , r coth r — 1 , , 

= , ■ ,2 (smh2r-2r), 
4sinh T 

Mr) = ^^f|^(-oshr + l), 



F(r) 



sinh T — r 
2 sinh r 

(2.16) 



In the far infrared, through the chain of Seiberg duahties that drop the size of the gauge 
group factors by M units, the D3-brane goes to zero and only the M fractional D3-branes 
remain. Thus, at the bottom of the cascade we are left with a pure A/" = 1 Yang-Mills 
theory with M isolated vacua. In this letter we will use the above supergravity description 
of this theory, proposed by Klebanov and Strassler, to study the glueball mass spectra of 
this theory by solving numerically the equations of motion describing the propagation of 
supergravity fields along the worldvolume of the branes. 

3 Supergravity Equations 

In this section we will write down the dimensional reduction of the linearized equations 
of motion of type IIB supergravity in the background described in section ^ The bosonic 
sector of the type IIB supergravity multiplet contains a graviton gp_c, a dilaton $, a 
zero-form R-R field C, two tensors -the NS-NS and R-R fields -B^i> and C/ip- and a 
R-R four-form C^^fi)- ■'■^ present notation carets denote indices that run over ten 
dimensions, greek indices denote four dimensional space and latin indices run over the 



internal space. Expanding the type IIB equations of motion p6| in background (dotted 
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fields) and fluctuations we get, 



Df^d^^ = -G-^,,G<'"\ (3.1) 
D'G,,, = -—[F,,^,,G''' + F^,^,,G''']+d^<!>G,,,, (3.2) 



- \g,^G'''G,-,,l (3.3) 
p = —s pi^a -mi f Q 4^ 

The presence of a non zero three-form in the background implies the coupling of the 
dilaton, metric and two form. In general this should require solving the system of cou- 
pled equations. However, these equations can be simplified by expanding the excitations 



in spherical harmonics p^l- Formally, one should expand in harmonics over the exact 
background. Nevertheless, since we are interested in calculating glueball masses, i.e., 
in confining effects, and these occur near the bottom of the cascade where the the S"^ 
shrinks to zero but the does not we can expand in Kaluza-Klein modes over the S^. 
The expansion for the dilaton is, 

$(x,y) = 5^0(x)'=r(y)^ (3.5) 

while for the two-form 

V(x,i/) = ^aJUa;)!''^!/), 
A,m{x,y) = J2[al%x)Yi^{y)+a';{x)D^Y'^{y)], 
AUx,y) = J2i^'^'(^^^n]iy)+^''i^)D[mY^fiy)], (3.6) 

where x = x'^ are four dimensional worldvolume coordinates are y = x^ are coordinates 
over the three sphere. The Lorentz type gauges D'^Amn = 0, D^Amp = can be used in 
order to fix a^^ = a^'^ = in ( ^.61) . 

First, we will consider the equation (3.1) for the dilaton which only requires scalar 
spherical harmonics. Substituting the expansion (|3.5|) into the field equation for the 
dilaton we see that the s-wave of the kinetic side becomes D^d^(l){x) Y{y), where we have 



omitted the /c = index. However, the three-form of the Klebanov-Strassler solution [|T^ 
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lives only on the internal coordinates so its spherical harmonic expansion involves Yj^n] 
but no scalar harmonics. Thus, the s-wave equation for the dilaton becomes D^d^(j){x) = 
or, explicitly, 

±d,[^d,<j>gn = 0. (3.7) 

In order to study the propagation of the supergravity two-form along the worldvolume 
of the branes we should consider the D°'Ga^u part of equation (3.2). Using 

c„ = € (^) w^.^/' ^y) (3-8) 

together with ( p.5| ) and ( |3.6| ) the s-wave equation of motion for the two-form in the 
worldvolume is 

^^^.M'W^aiv^^KVv"]^?""'^?^"'^'^?'^"^'] = -^^F,,,„,d'^<^a^^\ (3.9) 

where [ ] denotes antisymmetrization with strength one, and a^^y is the complexified two- 
form. 



4 Glueball Mass Spectra 

In this section we will study the discrete spectrum of glueball masses arising from propa- 
gation in ( p.8| ) of the type IIB dilaton field, and the complex antisymmetric field a^^. The 



SU{N + M) X SU{N) conifold theory has a U{1)r global symmetry |jT8|. In the present 
case, the deformation of the conifold breaks the U{1)r to Z2. Therefore, one expects 
a massless glueball in this theory. However, this massless glueball will couple to some 
combination of supergravity fields^, and thus we do not expect to find it in our analysis 
of the glueball spectra. 

In p|, ^ the correspondence between the AdSr, x background and primary chiral 
fields correlators was made explicit. This correspondence should, in principle, be modified 
for the deformed conifold background since the space is not asymptotically AdS anymore. 
However, as argued in |]TU||, only operators with A < |M can propagate all the way to 



r = 0. For this type of operators one expects that the correspondence should not be 



^We thank I. Klebanov for comments on this point. 
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greatly modified. Thus, one expects that the dilaton will couple to the dimension four 
operator Tr Ffj_^F'^'^, so that the spin zero glueball masses, which can be derived from the 
two point function of this operator, will arise from propagation of the dilaton. Similarly, 
the supergravity two-form should couple to a dimension six operator, i.e., to the spin one 
glueball. 

In order to obtain a pure glue gauge theory we should take the QgM limit. In this 
limit the curvature of the space is large and we are outside of the region where supergravity 
has small corrections. This situation is not new when dealing with string duals of pure 



glue theories. Thus, in the spirit of ||2T| , we proceed with the calculation and expect 
the corrections to masses to be small. With this caveat let us now consider the equation 
for the dilaton. Expanding the dilaton field in plane waves, 0(r, x) = /(r)e*'^'^, so that 
a mode of momentum k has a mass = —k"^ and using the metric (|2.8| ) the dilaton 
equation (|3.7|) becomes 



d 



dr 



df 



dr. 



3.21/3 " (sinh(2r) - 2Tf/^-f- - (fc^e^/^^) sinh2(r)/i(r)/ = 0. (4.1) 



As a boundary condition we will require that near the origin the function / must be 
smooth, so that df /dr = a.t t = 0. The asymptotics of /(r) as r — oo is obtained by 
demanding normalizability of the states. Since for large r, 

^ = 2-'3~'e^ sinh2(r)/i^/2^r) ~ re^/^r ^4 2) 

convergence of the integral / |$p signals an exponential behaviour of the solution near 
infinity, / ~ e""^, as can be suspected from direct inspection of (^4.1| ). Changing variables 
to /(r) = ■?/'(r)e"'^, the wave equation (|4.1|) becomes 

6[sinh(2r) - 2r]^/^^" + [6n2^/3(sinh(2r) - 2t)^^^ 

+ 2^/=^(cosh(2r) - l)(sinh(2r) - 2t)-'^/'^] f + [3n^2^/^{smh{2T) - 2r)2/3 (4.3) 
+ 2^/='n(cosh(2r) - l)(sinh(2r) - 2t)-^/^ - (kh^/^) sinh2(r)/i(r)] i/j = 0, 

This equation can be solved for large values of r by ipir) = Cie"""^ + C2e~^^^"'^^^^ which 
implies /(r) = ci + 026"^/^"^. Normalizability requires Ci = and we can fix C2 = 1. 
Thus, at infinity / ~ e~^^^'^. Equation ( [4.1D becomes then an eigenvalue problem; we 
want to find the values of k for which the equation has solutions matching the boundary 
conditions at and 00. We have done this numerically using a shooting technique ||2l|: 
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boundary conditions at infinity are taken as initial conditions and k is adjusted so that 
when numerically integrating ([4.1|) / is also smooth at r = 0. We should note that, 
numerically, infinity is taken as some large value of r such that t ^ k"^. On the other 
hand, since the boundary condition requires that / goes exponentially fast to zero at 
infinity, there is a numerical bound that does not allow us to push the numerical infinity 
as far as we want. This explains why we are restricted to obtain only the first two 
eigenvalues. The glueball masses thus obtained are shown in Table 1. They are measured 
in units of e^^^ which, as explained in [jl9|, sets the four dimensional mass scale of the 
field theory. 



State 


(Mass)2 


0++ 
0++* 


9.78 
33.17 



Table 1: Mass (squared) of the spin zero glueball and its first excited state obtained from 
supergravity. 



It is interesting to note that by varying the origin of integration one can see confining 
effects occurring in a small region close to the bottom of the cascade. Past this small 
region the spectrum becomes continuous signaling a conformal behaviour. 

Next we will study the equation of motion for the complexified two-form a^jy. Using 
again the metric ( p.8|) and the five- form fiux through the worldvolume, ( p.l3| ), equation 
r9D becomes 



3.2^/^ — 
dr 



/i(r)(sinh(2r) 



^ dr 



{kh^/^)smh\r)h\T)f] 



6lK 



15 



9sMH{T)e 



-2 



dr 



-e 



(4.4) 



where we have impossed solutions to be of the form a^^^r, x) = 



re 



ik.x 



and summation 



is implied over the a and (3 indices. We should note that as a consequence of the non 
constant five-form background of the solution of Klebanov and Strassler, the second order 
differential operator ( p.9|) can not be factorized into two first order operators, as in p7 . 



Thus, we should take care explicitly of the derivatives in the right hand side of ( [3.91 ). This 

(r) 



requires the complex decomposition ffj,u{T) = h^y{T) +ic^u{T). Impossing the gauge fixing 
condition a^,- = 0, so that the two- form has no components along the radial direction. 
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equation ( ^^ ) becomes a system of coupled equations. Similarly as in the dilaton case, 
the glueball masses can be obtained by numerically solving this system. The values are 
shown in Table 2. 



State 


(Mass)2 


1 — 
1—* 


14.05 
42.90 



Table 2: Mass (squared) of the spin one glueball and its first excited state obtained from 
super gravity. 



5 Discussion 

In this paper we have computed the dilaton and two-form excitations in the deformed 



conifold background recently constructed by Klebanov and Strassler [19|. Via the corre- 
spondence between supergravity and field theory, this corresponds to determining masses 
for glueballs in the dual effective field theory. Unlike the AdS^ x case, the Klebanov- 
Strassler background contains non trivial three-form and five-form fluxes. Far in the 
infrared the three-form flux prevents the three-cycle of the base from collapsing. We ob- 
tained the linearized type IIB equations by doing a Kaluza-Klein decomposition on this 
5*^. We solved the resulting eigenvalue problem by numerically integrating these equa- 
tions. This method is exact and allows us to find the first two excited states for the 
dilaton and two-form fluctuations. A WKB approximation could also be used to solve the 



equations. This approach was used in ||2^-||3^ to determine glueball masses for several 



finite temperature supergravity models. A comparison of WKB results with those of |21 
where a numerical method was used, shows that agreement increases for excited states. 
This is to be expected since the WKB approximation improves for large values of the 
masses. This might be at the root of the small discrepancy, when comparing ratios of 



eigenvalues, between our results for the spin zero glueball and those in |23], where the 
WKB approximation was used. Due to the fact that the equations involve extremely 
divergent functions our numerical method is not able to find higher excited states and we 
cannot compare with [^] for states where the WKB approximation is more reliable. It 



would be interesting to compare with lattice results but we are not aware of any lattice 
computation for supersymmetric glueballs. 
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The present computations can be generalized to higher spin glueballs. It would also be 
interesting to study mixed supergravity states; one would then expect to find a massless 
glueball as a consequence of the breaking of the U{1)r symmetry. 
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